3R-  8  0-0427 


m 


A  FINITE  ELEMENT  METHOD  FOR 
FIRST  ORDER  HYPERBOLIC  SYSTEMS 

by 

hUtc.ke.il  Ltu  kin  ^  ^ 

Department  of  Mathematics 
The  University  of  Michigan 
Ann  Arbor ,  MI  48109 


DTIC 

SELECTEI% 

JUN  2  4  1980 1  ■ 


DISTRIBUTION  STATEMENT  A 
Approved  for  public  release; 
Distribution  Unlimited 


AMS  (MOS)  subject  classification  (1970).  Primary  65N30,  65M15 

^  supported  by  AFOSR  under  Contract  F49620-79-C-0149  and 
by  a  Faculty  Research  Fellowship  from  the  Horace  Rackham 
School  of  Graduate  Studies,  The  University  of  Michigan. 


80 


9  14 1 


Abstract 


/ 

A  new  finite  element  method  is  proposed  for  the  numerical 
solution  of  a  class  of  initial-boundary  value  problems  for 
first  order  hyperbolic  systems  in  one  space  dimension.  An 
application  of  our  procedure  to  a  system  modeling  gas  flow 


in  a  pipe  is  discussed.  Asymptotic  error  estimates  are  derived 
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1.  Introduction. 

We  propose  and  analyze  a  new  finite  element  method  for 
the  numerical  solution  of  a  class  of  initial -boundary  value 
problems  for  first  order  hyperbolic  systems  in  one  space  dimen¬ 
sion.  Our  method  is  based  on  a  procedure . given  by  Platzman  [9]. 

A  generalization  of  our  procedure  for  problems  in  two  space 
dimensions  will  be  treated  in  a  later  paper  [8] . 

We  consider  problems  in  one  space  dimension  of  the  form 


■# 


%■ 


K 

c 


ut  +  (ai2(x,t)  v)x 


f1 (x,t,u,v) , 


(x,t)  e  [o,i]x [o  ,t] , 


(1.1) 


vt  +a2i(x't/U»v)ux  +  a22  (x,t,u,v)  vx 


*  f2(x,t,u,v),  (x,t)  e  [0,l]x[o,T], 


u(x,  0)  =  Ug ( x) , 


v(x,0)  =  Vg(x) ,  x  e  [0,1], 


v(0,t)  =  gQ(t) , 


v ( i , t )  =  gi(t) ,  t  e  [o,t] . 


We  assume  that  (u(x, t) , v(x, t) )  is  a  smooth  solution  of  (1.1) 


Let 
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l' 

k  . 

* 

:  r 

set 

{ (x , t ,  u ( x , t ) ,  v (x, t) )  I  (x,t)  e  [0,1] x [0,T] }  , 
and  assume  that  there  exists  a  positive  constant,  a  ,  such  that 


a12(x,t)  >  a, 


'•s 


(1.2) 


a21  (x,t,sirs2)  _>  a 


-2- 

for  all  (x,t,s^,s2)  e  CL  .  We  assume  also  that  a^  (x,t,s1#s2) 
and  f^(x,t,s^,s2)  are  Lipschitz  continuous  functions  of  their 
arguments  for  (x,t,s^,s2)  e  CL  • 

Initial-boundary  value  problems  of  the  form  (1.1)  occur, 
for  example,  by  scaling  the  space  variable,  x  ,  of  the  following 
first-order  system  modeling  the  transient  behavior  of  isothermal 
gas  flowing  in  a  pipe: 

pt  +  Gx  -  o,  (x,t)  e  [0,L]X(0,T], 

Gt  +  (a2  -  G2p“2)px  +  2Gp"1Gx 

(1.3)  *  -flGjGp'1,  (x , t )  e  [0,L]*(0,T]  , 

p (x, 0)  *  Pg(x),  G(x, 0)  »  Gq ( x ) ,  x  e  [0,L] , 

G(0, t)  *  gQ(t) ,  G(l » t)  *  g^ (t) ,  t  6  [0,T] , 

where  p  is  mass  density,  G  is  momentum  density  (averaged 
through  the  pipe  cross-section),  a  *  a(p)  is  the  isothermal 
speed  of  sound,  L  is  pipe  length,  and  f  =  f(|G|)  >0  is  a 
friction  factor.  We  assume  that  the  friction  factor  is  described 
by  the  Moody  diagram  (see  [12,  p.  288-289]).  In  this  case  there 
exist  positive  constants  Gc  and  fQ  such  that  f(|G|)  *  fQ|G|  1 
for  | G |  <  Gc  .  There  also  exists  a  positive  constant  f^ 
such  that 

lim  f  ( |  G I )  =  f.  . 

|G|—  1 

The  boundary  conditions  above  correspond  to  supplying  the  mass 
rate  of  flow  at  x  *  0,L.  Conditions  on  the  data  and  the  friction 


t 
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factor  guaranteeing  the  existence  of  global  smooth  solutions 
to  (1.3)  have  been  given  by  the  author  in  [7]. 

We  now  describe  the  finite  element  spaces  used  for  our 
procedure  for  (1.1).  Set  I  *  [0,1]  and  for  E  Cl  define 


Pk(E) 


{z:  I  •*>  ]R 


is  a  polynomial  of  degree  <  k}  . 


n 


i 


For  the  partition  6  -  {0  *  xQ 

<  X^  ...<  Xjj  =  1} 

we  define 

1i  =  ^i-l'*!1'  hi  “  xi“xi-l' 

and  h  =  max  h.  . 

Set 

(1.4)  ^k(r,6)  »  {z  e  Ck(I) 

z  e  *rUi) ,  i  *  1, 

2, . . . ,N}  . 

We  shall  often  write  '1/1  for 

^?k(r,6).  We  also 

set 

b'7*  *  ^^k(r'6)  * 

{||  |  z  e  ^?k(r,6)} 

• 

We  shall  assume  that  the  families  of  spaces  ^?k(r,6)  con¬ 
sidered  in  this  paper  are  based  on  meshes  6  that  are  quasi¬ 
uniform,  i.e.,  there  exists  >  0  independent  of  h  such 
that 

hi 

(1.5)  min  "JJ"  i  ci  • 

If  f ,g  e  L2 ( I) ,  denote 

[1 

<f,g>  *  fg  dx  . 

>  0 


We  propose  the  following  method  to  approximate  (1.1): 


-4- 


Pind  U:  [0,T]  -*■  7ft  ,  V:  [0,T]  ♦  ^  77L  such  that  (sup¬ 
pressing  explicit  dependence  on  x,t) 

<Ut»X>  -  <a12  V/Xx>  +  a12(l/t)g1(t)x(l) 

(1.6a) 

-  a12(0,t)g0(t)x(0)  -<f1(U,V)  ,x>,  X8^. 

<vt  +  a2l(U'V)Ux  +  a22(U'V)Vx' 

+  8[V(l,t)  -  gx(t)]  x(l) 

(1.6b)  +  8[V  (0 , t)  -  gQ(t) ]  X (0) 

-  <f2(u/v)/X>  / 

U(0)  *  uQ,  V(0)  S  v0  . 

The  terms  multiplied  by  the  constant  6  in  (1.6b)  repre¬ 
sent  penalty  terms  to  impose  the  boundary  conditions.  Let 

(1.7)  y  m  max  J  a22 (x,t,u(x, t) ,  v(x,t))|  . 

(x,t)e[0,l]x[0,T] 

We  show  that  the  scheme  (1.6)  is  convergent  if  8  >.  y/2.  If 
a22  5  that  the  scheme  (1.6)  is  convergent  for 

8  5  0.  We  must  also  require  that  >7!C  C^(I),  i.e.,  k  ^  1, 
if  a22  f  0.  Otherwise  the  term 

<a22(u,v)vx/  x>  ,  x  e  ^  , 


will  not  be  defined. 


We  denote  by  HJ ,  for  j  a  positive  integer,  the  Sobolev 

2 

space  of  functions  on  1  with  j  derivatives  in  L  (1)  and 


norm 


11*11?  -  J  f1  ^Tc  d*'  INI2 

J  le«0  JO  dx*  dx* 


k-0  JO  dx  dx 


z\\l  • 


We  shall  prove  the  following  theorem  in  Section  5: 


Theorem  1. 


Let  (u,v)  be  the  solution  to  (1.1)  and  assume 


that  there  exists  C2  <  ®  such  that 


sup  (||u||  +  II V||  )  <  C  , 

0<t<T 


{  [Klli  +  K«r]  i  C2  • 


Assume  that  U(0)  e  V(0)  e  ^satisfy,  for  some  C3  <  °°  , 
(1.8)  ||  U(0)  -  uQ  ||  +  ||  V(0)  -  vQ  1 1  <  C3hr  . 

Suppose  that  8  >_  y/2 ,  that  k  >_  1  if  a22  jL  0,  and  that  r  2 
if  the  system  (1.1)  is  nonlinear.  Then  there  exists  hQ  >  0  and 
C  such  that  the  solution  (U,V)  of  (1.6)  exists  on  [0,T]  for 
h  £  hg  and  such  that 


(1.9)  ||  U(t)  -u  (t)  ||  +  ||  V(t)  -v(  t)  ||  <  Chr  for  te  (0 ,  T)  ,  h<h 


We  note  that  the  estimate  (1.9)  is  of  optimal  order  for 
v.  It  has  not  yet  been  determined  whether  the  order  of  this 
estimate  for  the  approximation  of  u  can  be  improved  in 
general.  The  optimality  of  the  result  (1.9)  is  discussed  in 
more  detail  in  Section  5. 

Our  procedure  should  be  compared  to  the  standard  finite 
element  procedure  [2],  to  be  discussed  in  Section  2,  for  approxi¬ 
mating  the  solution  of  first-order  hyperbolic  systems.  Although 
high  convergence  rates  can  be  proven  for  both  the  standard  finite 
element  procedure  and  our  alternative  procedure,  there  are 
important  qualitative  differences  in  the  solutions  they  produce. 
It  is  shown  in  Section  2,  that  the  numerical  solution  produced 
by  the  standard  finite  element  procedure  has  dispersion  proper¬ 
ties  unlike  those  of  the  exact  solution  of  the  differential 
equations.  This  behavior  has  been  noted  by,  among  others, 
Hedstrom  [5]  and  Platzman  [9].  We  show  in  Section  2  that 

for  a  model  problem  which  is  a  linearization  of  (1.3)  our 
proposed  procedure  yields  a  solution  with  dispersion  properties 
similar  to  those  of  the  solution  of  the  differential  equations. 

In  Section  3  we  discuss  the  qualitative  nature  of  the 
solution  of  (1.3).  We  show  that  for  short  pipes  (L  small)  or 
rough  disturbances  the  dissipative  term  -f|G|Gp  1  is  rela¬ 
tively  unimportant  and  the  solution  to  (1.3)  is  approximated  by 
the  solution  to  a  wave  equation.  However,  for  long  pipes  and 
mild  disturbances  the  dissipative  term  -f|G|Gp  *  is  shown  to 
be  very  important  and  the  solution  to  (1.3)  is  approximated  by 
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the  solution  to  a  diffusion  equation.  Our  numerical  solution 
is  shown  to  be  close  to  that  of  the  solution  of  the  standard 
finite  element  method  for  the  wave  equation  or  the  diffusion 
equation  under  the  appropriate  conditions. 

Section  4  gives  a  result  which  indicates  that  accurate  long¬ 
time  integration  of  systems  such  as  (1.3)  is  possible  with  our 
procedure.  We  note  that  accurate  long-time  integration  of  (1.3) 
is  not  possible  with  the  standard  finite  element  method  due  to 
the  accumulation  of  round-off  and  truncation  error. 
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2*  Dispersion  analysis  of  an  example. 
Consider  the  hyperbolic  system 


ut  +  vx 

-  0 , 

(x,t) 

e  [0,1] 

*[0,T], 

vt  *  ux 

*  0, 

(x,t) 

e  [0,1] 

x [0,T] , 

(2. 

1) 

v(0,t)  ■ 

V(l,t) 

»  0,  t 

e  [o, t] 

9 

u(x,0)  - 

uQ(x), 

v(x,0)  -  vq (x) , 

x  e  [o,i]. 

The 

solution  to  ( 2 

.  1)  is 

easily  constructed 

through  Fourier 

analysis. 

If 

00 

uQ(x)  *  ^  an  cos  mrx. 

Vx)  = 

CD 

V  b  sin  mrx, 
j  n  ’ 

then 


u(x,t)  *  l  A  cos  ( nirt  +  0  )  cos  mrx  , 

0  n  n 

00 

v(x,t)  3  [  A  sinfrnrt  +  6  )  sin  mrx  , 
t  n  n 

where  AQ  =  aQ,  eQ  *  0,  Ar  cos  0n  =  aR,  An  sin  9n  *  bn  . 

Now  consider  two  approximation  procedures  for  (2.1).  Let 

hJ  *  {z  e  h1  |  z (0)  -  z(i)  *  o}  , 

'tn  *  ^?k(r,6)  n  hJ  , 

and  define  the  standard  finite  element  approximation  (2]  to  be 
the  functions  U  :  [0,T]  +7%,  V  :  [0,T]  -  77 ^  such  that 


1 
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<ut,w>  +  <vx,w>  -  o  ,  W  e  y?i' 

(2.2)  <vt,w>  +  <UX,W>  -  0,  w  e  77?  , 


U{0)  s  urt /  V(0)  *  vft  . 


It  is  instructive  to  study  an  explicit  solution  of  (2.2) 
in  the  case  k  *  0,  r  *  1,  and  6  is  a  uniform  partition. 

For  N  a  positive  integer,  set  h  *  i  and  *  ih, 
i  =  0,...,N.  Define  the  interpolation  operator  :  C(I)  •+■  7P7 
^0(1,6)  by  the  relations 

6^2{xi)  *  z(x^)  for  i  *  0, . . . ,N. 


then 


(2.3) 


U(0)  =  l  a  P'tcos  mrx)  , 
0  n 

N  _ 

V(0)  *  l  b  (y  (sin  mrx)  , 
1  n 


U(t)  =  l  An  cos( rnt  +  0n)  63( cos  mrx)  , 
N 

V(t)  =  V  A  sin(r  t  +  0  )  t>>(sin  mrx)  , 
j  n  n  n 


wnere  A.  =  a.,  A  cos  0„  =  a  ,  A  sin  0„  =  b  .  for 
Oun  nnn  nn 

n  *  1, . . . ,N,  and 
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Note  that  r„  =  0.  Thus,  we  see  that 

U(  t)  *  <^(cos  Nirx)  ,  V(t)  *  $(sin  Nttx)  =  0 

is  a  non-constant  steady  state  solution  to  (2.2).  Furthermore, 
we  see  from  Graph  #  1  that  the  most  spatially  oscillatory  com¬ 
ponents  of  the  solution  of  (2.2)  have  a  low  frequency  ir.  time 
even  though  the  spatially  oscillatory  components  of  the  solution 
of  (2.1)  have  high  frequency  in  time. 

We  can  also  see  this  phenomena  for  (2.2)  and  general  spaces 
(r,6)  as  follows.  Since 

dim  ^  7H°  +  1  <  dim  • 

there  exists  a  non-constant  t  e  7??  such  that 

(2.5)  <Xx'2>  =  -<x*x>  =  x  e  7??°  • 

Hence 

(2.6)  U(t)  =  2,  V(t)  =  0 

is  a  non-constant  steady  state  solution  of  (2.2).  We  note  by 
(2.5)  that  z  must  be  a  highly  spatially  oscillatory  function. 
Our  new  procedure  to  solve  (2.1)  is  the  following.  Find 
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u  :  [0,T]  *  771=  T7lk( r,«),  V  :  [0,T]  -  ^  77}  such  that 

<ut,w>  -  -  o,  we  771  • 

(2.7)  <Vt,W>  +  <UX,W>  =  0,  We^^?r 


U(0)  =  uQ,  V(0)  *  v0  . 


we  note  that  by  elimination  of  V  in  (2.7)  we  see  that  U 
satisfies 


<utt'W>  +  <ux,wx>  =  o,  w  e  771. 


Hence,  our  procedure  reduces  in  this  special  case  to  the  standard 
method  for  the  wave  equation. 

An  explicit  solution  to  (2.7)  can  be  constructed  as  follows. 
Let  be  the  kth  positive  Rayleigh-Ritz  approximate  eigen¬ 

value  with  eigenfunction  U.  e  7/1  for  the  problem 


(2.8) 


-u"  =  Au,  x  e  [0,11 , 


u’  (0)  *  u'  (1)  =  0. 


Thus,  0  <  A^  <  A2  <  • • •  <  Am  and 


Xk<^Uk#W^>  =  <uk  'W*>  '  W  e  77?  , 

where  M  *  dim  7^-1.  Let  cok  =  .  Also,  set  u>0  =  0 

and  UQ  =  1.  If 


IKJEJW- 1 


then 


M 

U(t)  s  l  \  cos(u)kt  +  8^)^  , 


V(t)  3  I  \  sin(o)kt  +  9k)  (-a)"1  Uk)x  , 


where  AQ  *  aQ,  0Q  -  0,  Ak  cos  8k  ■  ak  ,  Ak  sin  8k  *  bk  . 

It  is  well  known  [11,  p.223]  that  >  k  ir  for  k  *  1,...,M. 
Hence,  we  do  not  obtain  non-constant  steady  states  for  (2.7)  , 
and  states  which  are  spatially  oscillatory  have  high  frequency 
in  time. 

In  the  case  TPJ  ■  ^0(1,S)  where  5  is  a  uniform  mesh 

of  size  ,  N  a  positive  integer,  we  can  take 


and 

(2.9) 


Uk  =  cos  kirx) 


m  *  i ku h.  ,,  2  *  2 

u>k  *  2  sin(-^— )  (1  -  y  sm 


kirhj-1/2  h-l 


A  comparison  of  the  graph  of  u>k  and  rk  with  the  graph  of 
\>k  =*  kir  (see  Graph  #1)  shows  the  superior  dispersion  relation 
given  by  scheme  (2.7).  Intuitively,  the  graph  shows  thjat  only 
about  one-half  of  the  degrees  of  freedom  in  the  standard  method 
are  useful  in  approximating  the  solution  ,  whereas  for  the  pro¬ 
posed  method  all  of  the  degrees  of  freedom  are  useful. 


3.  Application  to  the  modeling  of  gas  flow  in  a  pipe. 

For  high  frequency  disturbances  solutions  of  equation 
(1.3)  behave  very  much  like  solutions  of  the  usual  second- 
order  wave  equation  and  for  low  frequency  disturbances  they 
behave  very  much  like  solutions  of  a  heat  equation.  In  this 
section  we  indicate  the  sense  in  which  this  is  true  for  a 
linearized  version  of  (1.3)  and  examine  the  behavior  of  our 
numerical  method  in  these  two  limiting  situations. 

If  we  linearize  (1.3)  about  constant  mass  density  p  >  0 
and  momentum  density  G  ■  0,  we  obtain  the  system  (after  scaling 
the  length) 

(3.1)  pt  +  L"1GX  *  0  ,  (x,t)  e  [0,1] x [0,T]  , 

(3.2)  Gt  +  L“lo2px  *  “fG, 

(3.3)  P (x, 0)  =  Pq ( x) ,  G(x, 0)  *  GQ(x),  x  e  [0,1], 

A  A 

for  positive  constants  L,  a  ,  and  f  where  L  is  pipe  length, 
o  =  a(p)  is  the  isothermal  speed  of  sound,  and  £  is  a  constant 
friction  factor.  We  consider  for  simplicity  the  case  of  homo¬ 
geneous  boundary  conditions 

(3.4)  G (0 , t )  «  G (1 , t )  *  0,  t  e  [0,T] . 

We  obtain  by  eliminating  the  variable  G  from  (3.1)-(3.2) 


the  damped  wave  equation 
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Ptt  -  l"232p  -  -bt  ,  (x,t)  e  [o,i] x [o,t] 


(3.5) 


p(x,0)  »  pQ (x) ,  p. (x,0)  »  -L  Gq  (x) ,  x  e  [0,1], 

U  t  X 


Px<0,t)  -  Px(l,t)  *  0, 


t  €  [0,T] . 


Now  suppose  that  the  initial  data  is  such  that  for  m  a 
positive  integer 


(3.5a) 


p(x,0)  =  a  cos(mirx),  p.(x,0)  =  8  cos (mux), 
m  z  in 


Then  it  is  easily  checked  that  the  solution  to  (3.5)  is  given  by 

(3.6)  p(x,t)  =  J~c*  e  m  +  c  e  m  |  cos (mux) 

|_  m  m  — * 


where 


± 

m 

f 

2 

L1  +  * 

± 

-y>, 

m 

Ym 

“  Yin 

h  -  4m2,2  If2  S2  f'2 


If  Y~  *  v+  for  some  m,  then  small  modifications  of  the 
1  m  'm 

following  arguments  are  necessary.  Let 

+  Y» 

p,  (x,t)  »  c  e  cos(mux), 

1  m 

- 

p-(x,t)  *  c  e  cos (mux). 

2  m 

Note  that  p  =  p^  +  p2  • 


m 
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We  first  consider  the  case  when 


(3.7) 


2  2  _-2a2  %-2 
ID  IT  L  O  f 


>>  1. 


A  simple  calculation  shows  that  under  the  condition  (3.7) 
we  have  that 


(3.8)  |fp.|«|L"Vpi  |, 

t  XX 

for  (x,t)  e  [0,l]x[0,T],  i  »  1,2. 
<p,<l/  e  C(  [0,1]  X  [0/T] )  we  say  |<j)|  << 
if  and  only  if  |4>|/|t|;|  <<  1  for 

it  follows  that  the  solution  to  (3. 
solution  to  the  wave  equation 


(For  functions 

1*1  for  (x,t)  e  [0,1] x [0/T] 
(x,t)  e  [0, 1] x [0,T] ) .  Hence, 
5)  is  approximated  by  the 


(3.9) 


0«  •  L  0  *  0 • 
px(0/t)  »  px(l,t)  *•  0, 

p (x, 0)  =  a  cos(nnrx)  , 
r  m 


(x,t)  e  [o , i] x [o ,t] , 
t  e  [o,T], 

p.  (X/0)  *  ft  cos(mirx)  ,  x  e  [0,1]. 
t  m 


Next,  we  consider  the  case  when 
(3.10)  m27r2  L~2d2  f“2  <<  1,  |6m|<<  |amf| 


Then  it  is  easily  verified  that 


(3.11)  Ip.  I  «  |l"292p  |,  I  ip  I  ,  for  (x ,  t )  e  [o,i]x[o,T]. 

2tt  2XX  2t 


It  also  follows  from  (3.10)  that 


so  that 


(3.12) 


for  (x,t)  e  [0,l]x[0,T]. 

Hence,  we  see  from  (3.11)  and  (3.12)  that  under  the  conditions 
(3.10)  we  have  that  the  solution  of  (3.5)  is  approximated  by 
the  solution  of  the  diffusion  equation 


“L“2a2pxx  =  -fpt  /  (x,t)  e  [0, 1] x [o,T] , 

(3.13)  px(0.t)  »  px(l,t)  »  0,  t  e  [0,T], 

p  (x, 0)  *  a  cos(mirx),  x  e  [0,1]. 

m 

We  shall  now  show  that  under  the  conditions  (3.7)  our 
numerical  solution  to  (3.1) -(3. 3)  is  close  to  the  standard 
finite  element  method  for  the  wave  equation  (3.9)  and  that  under 
the  conditions  (3.10)  our  numerical  solution  to  (3. 1) - (3. 3) is  close  to 
the  standard  finite  element  method  for  the  diffusion  equation 

(3.13) . 

Our  finite  element  solution  to  (3.1)-(3.3)  is  p:  [0,T] 

■*  such  that 


G;  [0,T] 
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(3.14b)  <Gt»x>  +  <l_1$2px*x>  -  x  e  -^7Pi  , 

(3.14c)  p (0)  *  pQ,  G(0)  =  GQ  . 

We  can  eliminate  G  from  (3.14)  to  obtain  the  following 
finite  element  equation  for  p  :  [0,TJ  ■+■ 

(3.15)  <Ptt'X>  +  <l_232px^  xx>  «  -<fpt*  x>  ,  x  e  in  . 
where  p(0)  is  as  in  (3.14c)  and  pfc(0)  e  7ZI  satisfies  by  (3.14a) 

(3.16)  <pt(0),  x>  *  <l“1g(0),  xx>  ,  x  e  yn  • 

As  in  our  discussion  of  the  differential  problem,  we  assume  that 


p<0)  *  amUm,  P4.(0)  -  8.0. 


m  m 


(we  use  here  the  notation  for  the  eigenproblem  (2.8)  introduced 
in  Chapter  2).  The  solution  to  (3.15)  is  then  given  by 


(3.17) 


where 


,  r*t  _  . _ 

»  ,  »  f  W  Ift  HI  |  .. 

p(t)  =■  t  Cm  e  +  cm  e  1  u, 


-  r_t 


m 


m 


-  |  [l  tJ  1  -  4xkL'V  f-2  ]  , 


®m 


+ 

m  m 


m 


r_  -  r 


m  m 


ai*  . . - 


»  4  •*  t 


J 
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An  analysis  similar  to  that  given  for  the  differential 
problem  then  establishes  that  if 


(3.18) 


A  IT2*2  f-2 


m 


>>  1/ 


we  have  that  the  solution  to  (3.15)  is  approximated  by  the 
solution  p:  [0,T]  -►  "7 71  to  be  standard  finite  element  approxi¬ 
mation  for  the  wave  equation 


(3.19) 


<ptt'x/>  +  a  2px»  xx^  *  °#  x  e  791 


Similarly,  if 


(3.20) 


\n  L-2  o2  r2  «  1,  |6m|  «  |c*mf | 


it  can  be  verified  that  the  solution  to  (3.15)  is  approximated 
by  the  solution  p  :  [0,T]  -*•  "T/i  of  the  standard  finite  element 
approximation  for  the  diffusion  equation 


^L~2  ®  2px'  Xx^  "  "f  ^pt'  x^  '  x  e 


(3.21) 


P(°)  =  a  U  . 
m  m 
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4.  Long-time  integration. 


We  shall  prove  the  following  result: 


Lemma  1.  Let  (p,G)  be  the  solution  to  the  continuous  problem 

A 

(3.1) -(3. 4) .  There  exists  positive  constants  C  and  D 
(depending  on  3,  L,  and  f ,  but  independent  of  T,  pQ,  GQ) 


such  that 


llp(t) 


(4.1) 


"  (  P0  +  II  G(t) 

JO 

<  =  «‘Dt(  II  «o  *  f  » 


°0  o*ll  *  II ®0H  1 


Also,  if  p:  (0 ,T]  -*■  772.,  G:  [0,T]  ia  the  numerical 

approximation  of  ( 3- 1) — (3- 4)  defined  by  (3.14),  then 


rl 

II P(t)  -  p (0)dx|j  +  |J G( t)  || 

1  0 

(4.2) 

<  Ce”Dt  (  ||  p  (0)  -  f1  p  (0)  dx||  +  ||  G(0)  ||  )  . 

Jo 


•  •  *  e  # 


Remark . 


This  result  implies  that  the  dependence  of  the  solution  p,  G 
of  (3-1)  —  (3.4)  at  tj^  >  tQ  on  pQ,  GQ,  and  the  boundary 
conditions  for  t  e  [0,tQ]  decays  exponentially  (except  for 

the  dependence  on  [  p(x,t«)dx,  of  course).  It  is  essential 

Jo  0 

that  a  numerical  procedure  for  (3.D  —  (3-4)  have  the  property 
(4.2)  if  the  accumulation  of  round-off  error  and  truncation 
error  is  to  be  prevented  from  destroying  the  accuracy  of  the 
numerical  solution  after  a  finite  interval  of  time. 

Proof  of  Lemma  1. 


We  shall  prove  the  estimate  (4.2)  for  the  numerical  solution 

given  by  (3.14).  The  proof  that  (4.1)  holds  for  the  differential 

problem  is  analogous.  Suppose  that  the  solution  p(t)  ■ 

M 

^  cm(t)Um  e  to  (3.14)  has  initial  conditions 

M  M 

p(0)  *  I  u  '  p.  (0)  -  j  s  0  . 

g  m  m  *t  g  wm  m 

We  assume  in  this  section  that  the  eigenfunctions  {Um>  have  been 

normalized.  Again,  if  y+  *  yZ  for  some  m,  then  small  modifications 

m  in 

of  the  following  arguments  are  necessary. 

It  follows  from  (3.17)  that 
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.  r*t  r“t 

Cm C t )  *  [o'**  e  m  +  c“e  m  J 
in  in  m 


By  the  orthonormality  of  the  set  {U  >  we  have  that 

in 


(4.3) 


f1  2  M 

||  p  (t)  -  p(0)dx||z  =  l  c  ( t) 

Jo  1  m 


It  follows  from  (3.14a)  that 


M 


(4.4) 


so 


(4.5) 


-1 


G(t)  =  Y  c'(t)  L A  *  U  , 
v  m  mm 

1  x 


M 


2  2-1 


I|g<«1!  -  S  =;(t> 


It  is  easily  checked  that  there  exists  a  positive  constant 
C  which  is  independent  of  m  and  5  such  that 


(4.6) 


i  +, 2  .  -, 2  2  .  „2  .2  -1. 

c  +  <=_  <  C  ( 0£_  +  S_  L  A  . 

1  m1  1  m 1  —  m  m  m 


We  also  note  that  since  >  A.  >  ir  for  all  partitions  6  and 

m  —  i  — 

for  m  =  1,...,M  it  follows  that 


(4.7)  Re  r'  <  Re  < 
m  —  1 


-|  [i  -7i  - 1,2#2  l-2  r2J 


f  =  D 


for  m  *  1, . . . ,M. 

Hence,  we  have  that 


(4.8) 


cm(t)  <  2e-2Dt  [|c;f2  +  |cj2] 


„  -2Dt  ,  2  .  „2  ,2  -1  , 
i  ce  1  “m  +  L  Am  1 
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Since  cm(0)  *  am  and  e^(0)  ■  8m  /  ifc  follows  from  (4.3)  and 
(4.5)  that 

(4'9)  l  “m  “  IIP  CO)  -jJ* (X0)<lac||2,  £  6*  L2  X^1  *  ||G(0)||  2  . 

Thus,  if  we  sum  (4.8)  for  m  =  1,...,M  we  obtain 


(4.10)  ||  p  (t) 


fl 

p  (0)dx|| 
JO 


2  <  C  e-2Dt 


(  II  p<0)  -  [1p(0)dx||  2  +  ||  GCO)  II2)  , 
JO 


Next,  it  is  easy  to  show  that  there  exists  a  positive  con¬ 
stant  C  which  is  independent  of  m  and  6  such  that 


(4.11) 


l'rm>2  Cl  i  C- 


Hence,  we  have  that 


||G(t)||2  -  I  c;<t)2  L2  a;1 


m  r  t  r_t  o  o  i 

=  1  (Cm  rm  e  +  Cmrm  e  >  2  L2  x'1 
m  m  mm  m 


(4.12) 


<  C  e 


-2Dt 


l  * 


+| 2  +  | c“ | 2 ] 
m1  1  m1  1 


<  C  e 


-2Dt 


(||  p  (0)  -  f  1p(0)dx  1 1  2  +  ||  G ( 0 )  1 1  2 )  , 
Jo 


We  note  that  the  estimate  (4.2)  is  not  valid  for  the 
standard  finite  element  solution  to  (3.1) -(3.4).  The  standard  fi¬ 
nite  element  solution  to < 3 . 1) - ( 3. 4 )  is  obtained  from  (2.2)  by 

A 

replacing  U  with  p  and  V  with  G  and  by  adding  -f  <G,W^ 
to  the  right  side  of  the  second  equation.  Since  (see  (2.5)) 
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..*.X  >  ■  *.  ' 


p(t)  =  z  ,  G(t)  S  0 

is  a  non-constant  steady  state  solution  it  is  clear  that  (4.2) 
cannot  be  verified  for  the  standard  finite  element  method. 


-24- 


5.  Proof  of  Theorem  1. 


We  shall  fix  r  >  k  ^  1  and  derive  optimal  order  error 
estimates  for  the  scheme  (1.6)  for  the  family  of  spaces 

7??  =  ( r ,  6 ) 


satisfying  the  quasi-uniformity  assumption  (1.5).  We  note  in 
the  statement  of  Theorem  1  that  we  may  allow  k  =  0  if  a22  =  0 
and  r  =  1  if  (1.1)  is  linear.  These  special  cases  can  be  proven 
by  a  variant  of  the  following  argument. 

In  what  folllows  we  denote  by  W],“  the  Sobolev  space  of 
functions  with  j  weak  derivatives  in  L°°(I)  and  norm 


'•'i  ■  Jo  IlSFlL-d)  ■  |z|  *  |2|°  • 

It  is  well-known  that  there  exists  C4  <  °°  ,  depending  only 
on  of  (1.5)  and  r,  such  that  the  following  inverse  hypo¬ 

thesis  holds  for  x  e  7%  > 


(5.1) 


»  '  Ixl  *  »  llxxll  <  c4  |l x  II  . 
h1/2  lx  *1  +  h  ||xxxll  <  =4  llxjl  • 


It  is  also  well-known  that  the  spaces  satisfy  the 
following  approximation  property: 

There  exists  C5  <  ®  such  that  for  2  <  s  <  r  +  1  and 
z  6  Hs  , 


tv. 
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(5.2) 


inf  (  ||  z  -  X  l|  +  h  ||  2  -  X  II  !  +  h  ||  2  “Xll2 

xe7 n 


*  h3/2|*-xli>  <  C5hS  H*lls  ' 


and  for  1  <  s  <  r  and  z  e  Hs  , 


(5.3) 


inf  ( ll2  “X  II  +  h||z -xll  +h1/2|z -xIq)  <C5h3  ||  z  H 

xe!^  1 


We  also  use  the  fact  that  there  exists  Cg  <  »  so  that  if  z  e  H3 
and  z  (0)  =  z  (1)  =  0  ,  then  we  can  choose  x  e  7?l  such 

X  X 

that  Xv(°*  =  Xvt1)  =  0  and 

X  X 


(5.4) 


2  -X||i  +  h  ||  X  ||  2  1  C6h  ||  z  ||  2 


In  what  follows,  C  will  denote  a  constant  which  depends 
on  the  Lipschitz  constants  for  a. .  and  f.,  a  »  and  C.  through 

^  J  X  X 

Cg,  but  which  is  independent  of  5  .  It  will  be  allowed  to  vary 
from  estimate  to  estimate.  When  the  arguments  for  a^ ^ ,  f^ 
are  omitted  we  assume  that  the  functions  are  evaluated  at 
(x,t,u(x,t),  v (x , t) ) . 

2 

We  wish  to  define  the  weighted  L  projection  of  v, 

R  *  R(v)  e  •  by  relations 

(5.5)  <a12(v  "  R(v))'W)  "  °'  w  e  7)\  . 

It  follows  by  the  approximation  property  (5.3)  and  (1.2)  that 
n 2  *  v  “  R(v)  satisfies 
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(5.6) 


n2  II  <  c  ||v||rhr 


It  then  follows  by  the  approximation  property  (5.3)  and  the 
inverse  hypothesis  (5.1)  that 


(5.7) 


hllnjllx  ♦  »1/2  ln2|  <  chr  ||v||r  . 


We  obtain  an  estimate  for  no  by  differentiating  (5.5).  We 

t 

then  obtain 


(5.8) 


^a12n2  '  +  dt  al2*  ^2'  = 


w  e  35 


it  follows  that 


(5.9) 


II n2  II  <  C  || n2 1|  +  c  inf  It vt  “  xll 

*  xe  ^7* 

<  Chr(  ||  V||  r  ♦  ||vt||r). 


As  before,  by  (5.3)  and  (5.1)  we  earn  then  obtain 

(5.10)  h||n2  ||x  +h1/2|r,2tl  <  chr(  ||v||r  +  llvtllr  >  • 

We  also  wish  to  define  the  following  approximation  of  u, 
Q  =  Q(u,v)  e  171  i  by  the  relations 

<a21(u  -  Q(u,v))x,xx>  +  <a22n2  '  Xx> 


(5.11) 


+  sn,(i,t)xv(i>  +  Sn,(0,t)xv(0)  *  o,  x  e 


J  (u  -  Q(u, 
Ja 


v) ) dx  =  0  . 


Set  nj_  *  u  -  Q(u,v).  It  follows  from  (5.1)  and  (5.7)  that 

|n2d.t)Xx(i)|  <  ch*-1  ||Xj[||  ||v||r,  x«»l. 

(5.12) 

|n2  (0,t)XxC0)  |  <.  chr’1||Xx||  ||v||r,  x  e  77L. 


Hence,  it  follows  that 

llnj.lt!  <  +Chr-1||v||r  +C  inf  |t  u  -  x  llj 

X  e77? 

(5.13) 

<  a>r_1  (  ||v||r  ♦  ||«||r). 

2 

We  can  derive  an  L  estimate  for  m  by  a  variant  of  a 

2 

frequently  used  duality  argument.  Let  f  e  1  (i),  f  dx  *  0 

JO 

and  determine  <p  by 


(a21^x^ x  =  f'  x  e  [0,11 ' 
»x<°>  =  <J>x(1)  *  °' 


Note  that 

♦x<*> 

Then  it  follows  that 


rx 

f  (s)ds. 

0 


IUII2  <  c  || « ||  . 


Also,  for  x  ey?t'  such  that  Xx(°)  =  Xx{1) 


0,  by  (5.11) 


JU  u  "m 


i 


4. 

;* 

•v, 

y 

£ 

i* 

$• 


* 
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<nx/  O  ^  nl,“  ta21*x*x^  "  ^a21T'lx'  *x> 


(5*14)  a  ^a21nl  7  x  ^  +  ^a22r,2  '  ^"X^x^ 

X  X 

+  <n2'  (a22<l>x)x>- 


Now  if  we  choose  x  to  approximate  $  as  in  (5.4)  so  that 

II*  “  X II  !  <  Cl>  II*  ll2  . 

it  follows  that 


(5.15)  |<n1.£>|  <  c  [h  Hnj.ll  x  +  h||  n2|[  x  *  ||n2|P  11*11  • 

Hence,  since 

,1 

n i  dx  -  o, 

Jo  *■ 


we  can  conclude  that 

(5.16)  Unjll  <  ar  t  ||u||r  +  ||v||r]  . 


We  can  derive  an  estimate  for  m  by  differentiating  (5.11) 

it 

with  respect  to  t  to  obtain  as  in  the  estimate  for  ||n1||^, 

(5-17)  IlnjJlj  <  Chr'l(  ||u||r  .  ||ot||r  +  ||v||c  *  ||vt||r). 

A  duality  argument  similar  to  the  previous  such  argument  establishes 


that 
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<5-l»  111!  II  <  Chr  (|| a|| r  +  ||v||r  +||aj|r  +||»t||r] 

t 


Set  ^  •  U-Q(u,v)  ,  52  *  V-R(v).  Then  from  (1.1),  (1.6a) 
and  (5.5)  we  see  that 


<^it'  x>  ‘  ^*12^2 *  xx  y 


(5.19) 


<nlfc.  X>  +  <Vu'v>  “  fi(u'v) »  x>'  x  e  »l' 


It  also  follows  from  (1.1),  (1.6b),  and  (5.11)  that 


<C2  +  »21(u'V)5l  +  a22(U'V)?2  'X> 

t  X  X 


(5.20) 


♦  SC2<1't)  X<D  +  B5 2<0,t)  x(0) 

■<Tn2  '  X^  +  <Ia21(u,v)  -  a21(u,v)lux/x^ 

t 

+  <[a22(u,v)  -  a22(0,V)  ]vx,  x  > 

+  <f2(u,v)  -  f  (u,v),  X>,  x  e  jrr • 


We  prove  the  theorem  by  a  variant  of  an  argument  used 
in  [2,6,10,13]  by  showing  that  there  exists  a  positive  constant 
C7  such  that  for  any  t  e  [o,T^  ,  if 


(5.21) 


(U,V)  exist*  on  [0,t] 


such  that 

{<x,s,  0(x,s),  V(x, s) )  |  (x,s)  e  [0,11  *  [0,t]>  £ 

and 

|  (u  -  0)X(S)  |  ♦  |  (V  -  V)x(s)  1  <  1  for  0  <  .  <  t. 

then 

(5.22)  ||  (u  -  0)  (s)  ||  +  ||  (v  -  V)  (s)  ||  <  C7hr,  0  <  s  <  t. 

Let  (B  be  a  compact  neighborhood  of 
{ (x, s ,  u(x,s) ,  v(x,s))  I  (X,s)  e  [0,1]  X  [0,T]} 
such  that  (8  interior  of  &  * 


It  follows  from  the  "inverse"  hypothesis  (5.1)  that  (5.22) 
implies  the  existence  of  a  positive  constant  C  so  that 

I  (u  -  U)  (s)  |  +  |  (V  -  V)  (s)  I  <Ch-  ,  o  <  s  <  t 

Hence,  there  exists  a  positive  constant  h^  such  that  if 
h  e  (0,  h1  ] ,  then  (5.2  2)  guarantees  that 

{ (x, S,  U(X,S),  V(X,S))  |  ( X , S )  e  [0,1]  X [0, t] }  <C  0  . 
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Similarly ,  it  follows  from  the  "inverse"  hypothesis  (5.1) 
that  (5.22)  implies  the  existence  of  a  positive  constant  C 
such  that 

r-| 

|  (u  -  U)x(s)  |  +  |  (v  -  V)x(s)  |  <  Ch  ,  0<s_<t. 

Hence,  since  r  >_  2  there  exists  a  positive  constant  h2  such 
that  if  he  (0,  hj  ]  ,  then  (5.22)  guarantees  that 

|  (u  -  U)x(s)  |  +  |  (v  -  V)x(s)  |  <  \  for  0  <  s  <  t  . 

Now  U(x,t)  ,  V(x,t) ,  |  (u  -  U)  (t)  |  ,  and  |(v-V)  (t)| 
are  continuous  functions  of  t.  Also,  by  the  assumption  (1.8) 
on  the  choice  of  U(x,0),  V(x,G)  and  the  above  discussion 
it  follows  that  for  h  _<  min(h^,h2)  we  have  (assuming  without 
loss  of  generality  that  C?  <  C3) 

{(x,0,  U(x,0)  ,  V(x,0))|  x  e  1}  O  S  , 

and 

I  (u  -  U)x(0)  |  +  |  (v  -  V)x(0)  |  <  j  . 

Suppose  we  have  proven  that  (5.21)  implies  (5.22).  We  can 
thus  conclude  from  the  existence  theory  for  ordinary  differ¬ 
ential  equations  that  if  h  £h0=min(  l^,  h2)  •  then  (U,V) 

exists  on  [0,T], 

{  (x,  s ,  U(x,s),  V  (x ,  s ))  |  ( x ,  s )  e  [0,1]  X  [0,T]  }  C  (B  , 

|  (u  -  U)x(s)  I  +  I  (v  -  V)x(s)  I  <  i  ,  0  <  s  <  T  , 
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||  (u  -  U)  (s)  II  +  II  (V  -  V)  (s)  II  <  C?hr,  0  <  s  <  T. 


So,  we  assume  that  (5.21)  is  valid.  We 
shall  proceed  to  prove  the  estimate  (5.22).  Since 


{ (x , s ,  U(x,s),  V(x,s) )  |<x,s)  e  (0,1]  x  [0,t]  C.  CL 


we  can  assume  (1.2)  and  the  Lipschitz  bounds  on  the  coefficients 
where  needed.  We  also  note  that  it  follows  from  (5.21)  and 
Sobolev's  inequality  that 


|ux(s)|  +  |vx(s)|  <  |ux(s)|  +  |vx(s)|  +  1 


(5.23) 


for  0  <  s  <  t. 


<  C  {  ||  u (s)  ||  _  +  ||  v ( s)  1 1  2  }  +  1  <  c  C2  +  1 


,s‘ 


We  set  ag.  e  77?  to  be  the  L*  projection  of  g. 

i  *12  J- 

onto  77l ,  i.e. , 


(5.24) 


^  ITT  '  X  e  • 


We  use  the  following  variant  of  a  lemma  of  Dupont  and  Wahlbin 

— *21 

[4]  to  estimate  the  error  ag.  -  -=i  r.  . 

i  a12  1 

Lemma  2.  There  is  a  constant  C  such  that  if  b  is  a  continuously 

2 

differentiable  function  on  [0,1] ,  $  e  7y,  and  y  is  the  L 

projection  of  b$  onto  »  then 

ll<b*  -  *>xll  <  C|bx|  |U||  . 

Proof  of  Lemma.  For  x  6  >*  '  we  have  by  (5-1) 

II  «  0>xl|  <  ||  <b*  -  X)J|  ♦  ||  <x  -  0>x|| 


(5.25) 


<  II  <b0  -  X>x||  +  Ch  ||  X  “  Oil 


<  l|(bO“X)J|  +  Ch“A(||  b*  -  xll  +  || b0  -*| 


It  follows  directly  from  [4]  that 


IlN  -*||  <  Ch  |b  I  |U||  . 


The  argument  in  [4]  also  shows  that  one  can  construct  x  6 


such  that 


||b*  -x||  +  h  ||  (b^  -  x>*ll  <  chibxl  IUII 


The  result  follows  by  substituting  the  above  two  estimates  in 
(5.25).  Q.E.D. 


We  apply  the  lemma  to  obtain 


(5.26) 


<a«l>x  -  iC  |<S^>x|  lUlll 


We  note  that  from  (1.2)  it  follows  that  there  exists  a  constant 
C  >  1  such  that 


c'1  lUlll  <  II  Sill  i  c  lUil 


We  now  take  x  =  a$i  e  77?  as  our  test  function  in  (5.19) 
We  obtain  since  e  y?l  that 
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(5.27) 


<V  hT  ^  ‘  *  ?1'  ^ 

_  l/c  i±  hi)  r  > 

2^1'  1  dt  a,  J  n' 


and  by  (5.26) 


(5.28) 


<ai2C2'  *  ^a12C2'  (a^ 


+  0(  U  1| 2  +  ll  11  2  ) 

=  <^a21C2'5l  >  +  °(  llC2[l2  +  (l^ll2) 


Hence,  we  can  obtain  from  (5.19),  (5.27)  and  (5.6),  (5.16), 


and  (5.18)  that 


(5.29) 


TdE<^i'  I^i>  "  <fa2i?2^ix> 


<  C  [  i|  Cxl| 2  +  ||C2I|2  1  +  Ch2r 


If  we  set  x  *  £2  e  3x  '7%.  in  (5*20)  and  use  (5*6)»  (5.9),  (5,16) 


and  (5.23)  we  obtain 


(5.30) 


'^2^  +  <sa21?l  '^2/^  +  <^a22^2  '^2^ 
t  x  x 

+  8C2(1»S)2  +  8£2 (Or  S) 2 


<  ch2r  +  C  [  lUi  il 2  +  ||  52 II 2  1 


~~r 


However 


^a22S2x/S2')  2  jQ  (a22?2,xdX"  2  <'(a22,x 

(5.31) 

<  |  Y(C2(l,s)2  +  C2(0,s)2)  +  C  || C2 II2  • 

Thus,  if  we  add  (5.29)  to  (5.30)  and  use  (5.31)  we  obtain 

i&Kh'  IT1  Sl>+  lic2||2  ]  +  (6  -JY)  (Z2a.s)2  +  c2(0,s)2) 

X  » 

(5.32) 

<  c  [  jiqii2  +  |U2||2  ]  +  ch2r  . 

We  now  obtain  from  Gronwall ' s  lemma  the  result 

(5.33)  ||S1(s>||2  +  IU2(s)  i|2  +  (B  -  |y)  *  U2  (1  .a) 2  +  £2  (0  ,a)  2)  da 

2r 

£  Ch  for  0  £  s  £  t. 

Since  it  follows  from  (5.6)  and  (5.16)  that 


Hii2  +  ll»2ii2  <  c»2t  , 

we  have  shown  that 


II  (U-U)  (S)  ||  +  ||  (V-V)  (S)  ||  £  Chr  f  0  £  S  £  t. 

Thus,  we  have  shown  that  (5.21)  implies  (5.22).  Hence,  we  can 
conclude  (as  our  theorem  asserts)  that  there  exists  h^  >  0 
such  that  if  0  <  h  <  h„  ,  then  (U,V)  exist  on  [0,T1  and 
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(5.34)  ||  (u-U)  (t)  ||  +  ||(v-V)  (t) ||  <  Chr  for  0  <  t  <  T. 

Q.E.D. 

We  also  note  that  the  following  estimate  is  valid  [1] : 

(5.35)  | n2 |  <  Chr  |v|r  . 

We  can  use  this  estimate  to  conclude  from  (5.33)  the  following 
corollary. 

Corollary  1.  If  in  addition  to  the  conditions  of  Theorem  1 
we  have  that 

T 

(5.36)  f  |v| 2  dt  <  a, 

Jo  r 

and  0  >  ^  ,  then  the  estimate 
T 

(5.37)  f  [  |  v  ( 1 ,  t )  -V(l,t)|2  +  | v(0, t)  -  V(0,t) |2]dt  <  Ch2r 

Jo 

is  valid. 

We  now  consider  the  optimality  of  the  estimate  (1.9). 

First,  we  consider  the  problem 

ut  +  vx=0,  (x,t)  e  [0,1]  x  [0 ,  T]  , 

vt  +  ux  -  -tv  , 

(5.38) 

v(0, t)  *  gQ(t)  ,  v(l,t)  *  gL(t),  t  e  [0,T]  , 

u(x, 0)  *  Uq (x) ,  v(x,0)  *  vo(x)f  x  €  t0,l]  , 
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where  £  is  a  constant.  The  next  lemma  shows  that  if  the 
solution  (u,v)  to  (5.38)  is  smooth  enough  and  if 

1|  uQ  -  U  (0)  ||  <  Chr+1, 

(5.39) 

<v0  -  V(0),  x >  =0,  x  e  ^7 n  . 

then  the  improved  estimate 

(5.40)  h_1  ||  u ( t )  -  U(t)  ||  +  ||v(t)  -  V(t)  ||  <  Chr,  te  [0,T]  , 

can  be  obtained.  We  note  that  the  estimate  (5.40)  is  of  optimal 
order  in  h  for  both  u  and  v. 

Lemma  3 .  Let  (u,v)  be  the  solution  to  (5.38)  and  assume 
that  there  exists  Cg  such  that 


Let  (U,  V)  be  the  solution  to  (1.6)  with  a12  =1,  =  °» 

A 

a21  =  1,  a22  =  0,  f2(u,v)  =  -fv  ,  8  =  0,  k  0,  and  r  >  1. 

Assume  that  (U(Q),  V(0))  satisfies  (5.39).  Then  the  estimate 

(5.40)  is  valid. 

Proof.  It  is  easily  checked  that (using  the  notation  of  Theorem  1) 
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(5.41)  ,  x>  -  <52,  xx>“<n1*x>*  x  e  771' 

t  t 

and 

(5.42)  <"52  '  X>  +  <5l  •  X>-  -K?2-X>-  X  e  ^  • 

We  set  X  *  in  (5.41),  x  *  C2  ^5*42^'  and  add  the 

resulting  equations  to  obtain 

i  as  i  IU2.112  +  iu2il2  J 

(5.43) 

-  7  II ni  II 2  +  7 II ^ill 2  +  Ul  I|C2II2  • 

Now  it  follows  from  standard  estimates  that  in  this  case 
tlnj.II  <  Chr+l  |M|rtl.  IJnjJI  <  Chr+1  ||ut||r+1  , 

(5.44) 

II  n2 1|  <  chr  ||  v  ||r  . 

Also,  we  can  conclude  from  (5.39)  and  (5.44)  that 

(5.45)  ||£1<0:||  <  Chr+1  ,  C2(0)  SO. 

Hence,  it  follows  from  (5.43),  (5.44),  and  (5.45)  and  Gronwall’ 
inequality  that 

(5.46)  ||  (t)  || 2  +  |U2(t)||2  <  Ch2(r+1),  t  €  [0  ,T]  . 

The  result  (5.40)  now  follows  from  (5.44)  and  (5.46).  Q.E.D. 


We  conjecture  that  the  estimate  (5.40)  is  valid  in 
general  if  the  solution  (u,v)  to  (1.1)  is  smooth  enough, 
if  a22  =  0,  and  if 

||u0-D(0)||  <  Chr+1, 

<a12<v0  '  V(0))'  X>  ’  °'  *  • 

Next,  consider  the  system 

ufc  +  vx  *  f1(x),  (x,t)  €  tO,l]x[0,T), 

Vt  +  Ux  +  Vx  "  f2(x)' 

(5.47) 

v(0,t)  *  g0(t)  ,  v ( 1 , t )  *  g^t),  t  e  [0,T]  , 
u(x, 0)  *  uQ(x),  v(x,0)  *  vQ(x),  x  e  [0/1]. 

Let  the  space  r/6)  where  6  is  a  uniform  partition 

with  mesh  length  h,  i.e.,  h  *  ^  for  N  a  positive  integer 
and  Xj  *  jh. 

Set  m  58  u  “  Q(u,v)  as  before.  Let  (U,V)  be  the 
solution  of  scheme  (1.6)  for  (5.47)  with  6*1  and  initial 
conditions 

(5.48)  V(0)  *  R(vq) ,  U(0)  *  Q(u0,vQ)  . 

Let  *  7?f^(2,5).  It  can  be  shown  that  if  (u,v)  is  a 

smooth  solution  to  (5.47),  then  there  exists  a  constant  C9 


such  that 
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(5.49)  |( nx ((  +  lln^l  £  Cgh5/2. 

This  result  can  be  used  in  the  analysis  of  Theorem  1  to  show 
that  if  (u,v)  is  a  smooth  solution  to  (5.47),  then  the  result 

(5.50)  h"1/2  || u (t)  -  U(t)  ||  +  || v (t)  -  V(t)  ||  <  Ch2  , 

t  e  ro,T]  , 

can  be  proven.  The  estimate  (5.50)  cannot  be  improved  since  it  can  be 

proven  by  techniques  similar  to  those  used  in  [3]  that  the  estimate  at 
t  »  0, 

h_1/2  ||  uQ  -  Q(u0,vq)  ||  +  ||v0-R(v0)||  <  Ch2  , 

cannot  be  improved. 

It  can  also  be  shown  that  if  •m-vn2  (4,6)  and  if  the 
initial  conditions  are  determined  by  (5.48),  then  the  estimate 
(1.9)  cannot  be  improved. 


t 
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